Modeling biophysical feedbacks in the Earth system to investigate a fire-controlled hysteresis of tropical forests

Tropische Regenwälder sind durch anthropogene Aktivitäten gefährdet und wurden als Kippelement identifiziert. Ein Kippen in einen neuen Zustand könnte tiefgreifende Auswirkungen auf das globale Klima haben, sobald die Vegetation von einem bewaldeten in einen Savannen- oder Graslandzustand übergegangen ist. Waldbrände können die Grenze zwischen Savanne und Wald verschieben und somit das dynamische Gleichgewicht zwischen diesen beiden möglichen Vegetationszuständen unter sich änderndem Klima beeinträchtigen. In der vorliegenden Doktorarbeit wurde ein neues Erdsystemmodell entwickelt und angewendet, um explizit die Auswirkungen von Feuer, Klimawandel und Landnutzung auf eine potenzielle tropische Hysterese abzuschätzen. In den ersten beiden Teilen der Arbeit wurde das Vegetationsmodell LPJmL vor allem in Hinblick auf Feuersimulation verbessert und anschließend biophysikalisch an das Erdystemmodell CM2Mc gekoppelt. Im dritten Teil dieser Arbeit wurde das resultierende Modell CM2Mc-LPJmL schließlich angewendet, um wichtige biophysikalische Feuer-Vegetations-Klima-Rückkopplungen und einen potentiellen Kipppunkt bzw. eine Hysterese der tropischen Wälder zu untersuchen. Die Ergebnisse der Experimente zeigten, dass eine alleinige Klima Störung nicht zu einem großflächigen Kipppunkt tropischer Wälder führt. Andererseits führte die vollständige Entwaldung bei einer erhöhten CO2-Konzentration von über 450 ppm und die Wirkung von Waldbränden zu einer Verschiebung großer Teile des Amazonas Regenwaldes in einen stabilen Graslandzustand. Die Leistung dieser Arbeit ist die Entwicklung eines neuen Erdsystemmodells, das die Vorteile des umfassenden dynamischen Vegetationsmodells LPJmL und eines prozessbasierten Feuermodells mit dem geringen Rechenaufwand von CM2Mc verbindet. Diese Doktorarbeit untersuchte zum ersten Mal den expliziten Einfluss von Feuer auf tropische Kipppunkte und auf eine mögliche vegetative Erholung in einem umfassenden feuerfähigen Erdsystemmodell.Tropical rain forests are endangered by anthropogenic activities and are recognized as one of the terrestrial tipping elements. An ecosystem regime change to a new state could have profound impacts on the global climate, once the biome has transitioned from a forest into a savanna or grassland state. Fire could potentially shift the savanna-forest boundary and hence impact the dynamical equilibrium between these two possible vegetation states under a changing climate. In this thesis, a new Earth system model was developed and applied to explicitly estimate the impact of fire, climate change and land-use on a potential tropical tipping point and hysteresis. The first part of this thesis describes the improvement of simulating fire within the dynamic global vegetation model (DGVM) LPJmL (Lund-Potsdam-Jena-managed-Land). In the second part, the improved LPJmL model was biophysically coupled to the Earth system model CM2Mc, which involved numerous changes in the original LPJmL model. In the third part of this thesis, the resulting model CM2Mc-LPJmL was finally applied to investigate important biophysical fire-vegetation-climate feedbacks and a potential tipping point and hysteresis of tropical forests. The results of the modeling experiments indicated that a sole climate disturbance does not lead to a large-scale tipping of tropical forests into a savanna or grassland state. On the other hand, complete deforestation alongside elevated CO2 above 450 ppm and the impact of fire led to a shift of large parts of the Amazon into a stable grassland state. The contribution of this thesis is the development of a new Earth system model, including the advantages of the comprehensive dynamic vegetation model LPJmL, a process-based fire model and the low computation cost of CM2Mc. This thesis studied for the first time the explicit impact of fire on tropical tipping points and a possible vegetation recovery in a comprehensive fire-enabled Earth system model

Physically Constrained Generative Adversarial Networks for Improving Precipitation Fields from Earth System Models

Precipitation results from complex processes across many scales, making its accurate simulation in Earth system models (ESMs) challenging. Existing post-processing methods can improve ESM simulations locally, but cannot correct errors in modelled spatial patterns. Here we propose a framework based on physically constrained generative adversarial networks (GANs) to improve local distributions and spatial structure simultaneously. We apply our approach to the computationally efficient ESM CM2Mc-LPJmL. Our method outperforms existing ones in correcting local distributions, and leads to strongly improved spatial patterns especially regarding the intermittency of daily precipitation. Notably, a double-peaked Intertropical Convergence Zone, a common problem in ESMs, is removed. Enforcing a physical constraint to preserve global precipitation sums, the GAN can generalize to future climate scenarios unseen during training. Feature attribution shows that the GAN identifies regions where the ESM exhibits strong biases. Our method constitutes a general framework for correcting ESM variables and enables realistic simulations at a fraction of the computational costs

Constraining modelled global vegetation dynamics and carbon turnover using multiple satellite observations

The response of land ecosystems to future climate change is among the largest unknowns in the global climate-carbon cycle feedback. This uncertainty originates from how dynamic global vegetation models (DGVMs) simulate climate impacts on changes in vegetation distribution, productivity, biomass allocation, and carbon turnover. The present-day availability of a multitude of satellite observations can potentially help to constrain DGVM simulations within model-data integration frameworks. Here, we use satellite-derived datasets of the fraction of absorbed photosynthetic active radiation (FAPAR), sun-induced fluorescence (SIF), above-ground biomass of trees (AGB), land cover, and burned area to constrain parameters for phenology, productivity, and vegetation dynamics in the LPJmL4 DGVM. Both the prior and the optimized model accurately reproduce present-day estimates of the land carbon cycle and of temporal dynamics in FAPAR, SIF and gross primary production. However, the optimized model reproduces better the observed spatial patterns of biomass, tree cover, and regional forest carbon turnover. Using a machine learning approach, we found that remaining errors in simulated forest carbon turnover can be explained with bioclimatic variables. This demonstrates the need to improve model formulations for climate effects on vegetation turnover and mortality despite the apparent successful constraint of simulated vegetation dynamics with multiple satellite observations

Improving the LPJmL4-SPITFIRE vegetation–fire model for South America using satellite data

Vegetation fires influence global vegetation distribution, ecosystem functioning, and global carbon cycling. Specifically in South America, changes in fire occurrence together with land-use change accelerate ecosystem fragmentation and increase the vulnerability of tropical forests and savannas to climate change. Dynamic global vegetation models (DGVMs) are valuable tools to estimate the effects of fire on ecosystem functioning and carbon cycling under future climate changes. However, most fire-enabled DGVMs have problems in capturing the magnitude, spatial patterns, and temporal dynamics of burned area as observed by satellites. As fire is controlled by the interplay of weather conditions, vegetation properties, and human activities, fire modules in DGVMs can be improved in various aspects. In this study we focus on improving the controls of climate and hence fuel moisture content on fire danger in the LPJmL4-SPITFIRE DGVM in South America, especially for the Brazilian fire-prone biomes of Caatinga and Cerrado. We therefore test two alternative model formulations (standard Nesterov Index and a newly implemented water vapor pressure deficit) for climate effects on fire danger within a formal model–data integration setup where we estimate model parameters against satellite datasets of burned area (GFED4) and aboveground biomass of trees. Our results show that the optimized model improves the representation of spatial patterns and the seasonal to interannual dynamics of burned area especially in the Cerrado and Caatinga regions. In addition, the model improves the simulation of aboveground biomass and the spatial distribution of plant functional types (PFTs). We obtained the best results by using the water vapor pressure deficit (VPD) for the calculation of fire danger. The VPD includes, in comparison to the Nesterov Index, a representation of the air humidity and the vegetation density. This work shows the successful application of a systematic model–data integration setup, as well as the integration of a new fire danger formulation, in order to optimize a process-based fire-enabled DGVM. It further highlights the potential of this approach to achieve a new level of accuracy in comprehensive global fire modeling and prediction

Climate-induced hysteresis of the tropical forest in a fire-enabled Earth system model

Tropical rainforests are recognized as one of the terrestrial tipping elements which could have profound impacts on the global climate, once their vegetation has transitioned into savanna or grassland states. While several studies investigated the savannization of, e.g., the Amazon rainforest, few studies considered the influence of fire. Fire is expected to potentially shift the savanna-forest boundary and hence impact the dynamical equilibrium between these two possible vegetation states under changing climate. To investigate the climate-induced hysteresis in pan-tropical forests and the impact of fire under future climate conditions, we employed the Earth system model CM2Mc, which is biophysically coupled to the fire-enabled state-of-the-art dynamic global vegetation model LPJmL. We conducted several simulation experiments where atmospheric CO2 concentrations increased (impact phase) and decreased from the new state (recovery phase), each with and without enabling wildfires. We find a hysteresis of the biomass and vegetation cover in tropical forest systems, with a strong regional heterogeneity. After biomass loss along increasing atmospheric CO2 concentrations and accompanied mean surface temperature increase of about 4 ∘C (impact phase), the system does not recover completely into its original state on its return path, even though atmospheric CO2 concentrations return to their original state. While not detecting large-scale tipping points, our results show a climate-induced hysteresis in tropical forest and lagged responses in forest recovery after the climate has returned to its original state. Wildfires slightly widen the climate-induced hysteresis in tropical forests and lead to a lagged response in forest recovery by ca. 30 years

CM2Mc-LPJmL v1.0: biophysical coupling of a process-based dynamic vegetation model with managed land to a general circulation model

The terrestrial biosphere is exposed to land-use and climate change, which not only affects vegetation dynamics but also changes land–atmosphere feedbacks. Specifically, changes in land cover affect biophysical feedbacks of water and energy, thereby contributing to climate change. In this study, we couple the well-established and comprehensively validated dynamic global vegetation model LPJmL5 (Lund–Potsdam–Jena managed Land) to the coupled climate model CM2Mc, the latter of which is based on the atmosphere model AM2 and the ocean model MOM5 (Modular Ocean Model 5), and name it CM2Mc-LPJmL. In CM2Mc, we replace the simple land-surface model LaD (Land Dynamics; where vegetation is static and prescribed) with LPJmL5, and we fully couple the water and energy cycles using the Geophysical Fluid Dynamics Laboratory (GFDL) Flexible Modeling System (FMS). Several improvements to LPJmL5 were implemented to allow a fully functional biophysical coupling. These include a sub-daily cycle for calculating energy and water fluxes, conductance of the soil evaporation and plant interception, canopy-layer humidity, and the surface energy balance in order to calculate the surface and canopy-layer temperature within LPJmL5. Exchanging LaD with LPJmL5 and, therefore, switching from a static and prescribed vegetation to a dynamic vegetation allows us to model important biospheric processes, including fire, mortality, permafrost, hydrological cycling and the impacts of managed land (crop growth and irrigation). Our results show that CM2Mc-LPJmL has similar temperature and precipitation biases to the original CM2Mc model with LaD. The performance of LPJmL5 in the coupled system compared to Earth observation data and to LPJmL offline simulation results is within acceptable error margins. The historical global mean temperature evolution of our model setup is within the range of CMIP5 (Coupled Model Intercomparison Project Phase 5) models. The comparison of model runs with and without land-use change shows a partially warmer and drier climate state across the global land surface. CM2Mc-LPJmL opens new opportunities to investigate important biophysical vegetation–climate feedbacks with a state-of-the-art and process-based dynamic vegetation model

Variable tree rooting strategies are key for modelling the distribution, productivity and evapotranspiration of tropical evergreen forests

A variety of modelling studies have suggested tree rooting depth as a key variable to explain evapotranspiration rates, productivity and the geographical distribution of evergreen forests in tropical South America. However, none of those studies have acknowledged resource investment, timing and physical constraints of tree rooting depth within a competitive environment, undermining the ecological realism of their results. Here, we present an approach of implementing variable rooting strategies and dynamic root growth into the LPJmL4.0 (Lund-Potsdam-Jena managed Land) dynamic global vegetation model (DGVM) and apply it to tropical and sub-tropical South America under contemporary climate conditions. We show how competing rooting strategies which underlie the trade-off between above- and below-ground carbon investment lead to more realistic simulation of intra-annual productivity and evapotranspiration and consequently of forest cover and spatial biomass distribution. We find that climate and soil depth determine a spatially heterogeneous pattern of mean rooting depth and below-ground biomass across the study region. Our findings support the hypothesis that the ability of evergreen trees to adjust their rooting systems to seasonally dry climates is crucial to explaining the current dominance, productivity and evapotranspiration of evergreen forests in tropical South America.Fil: Sakschewski, Boris. Potsdam Institute for Climate Impact Research; AlemaniaFil: Von Bloh, Werner. Humboldt-Universität zu Berlin; AlemaniaFil: Drüke, Markus. Humboldt-Universität zu Berlin; AlemaniaFil: Sörensson, Anna. Consejo Nacional de Investigaciones Científicas y Técnicas. Oficina de Coordinación Administrativa Ciudad Universitaria. Centro de Investigaciones del Mar y la Atmósfera. Universidad de Buenos Aires. Facultad de Ciencias Exactas y Naturales. Centro de Investigaciones del Mar y la Atmósfera; ArgentinaFil: Ruscica, Romina. Consejo Nacional de Investigaciones Científicas y Técnicas. Oficina de Coordinación Administrativa Ciudad Universitaria. Centro de Investigaciones del Mar y la Atmósfera. Universidad de Buenos Aires. Facultad de Ciencias Exactas y Naturales. Centro de Investigaciones del Mar y la Atmósfera; ArgentinaFil: Langerwisch, Fanny. Universitat Potsdam; AlemaniaFil: Billing, Maik. Universidade Federal de Santa Catarina; BrasilFil: Bereswill, Sarah. Universidade Estadual de Campinas; BrasilFil: Hirota, Marina. Potsdam Institute for Climate Impact Research; AlemaniaFil: Oliveira, Rafael Silva. Potsdam Institute for Climate Impact Research; AlemaniaFil: Heinke, Jens. Potsdam Institute for Climate Impact Research; AlemaniaFil: Thonicke, Kirsten. Potsdam Institute for Climate Impact Research; Alemani

Variable tree rooting strategies are key for modelling the distribution, productivity and evapotranspiration of tropical evergreen forests

A variety of modelling studies have suggested tree rooting depth as a key variable to explain evapotranspiration rates, productivity and the geographical distribution of evergreen forests in tropical South America. However, none of those studies have acknowledged resource investment, timing and physical constraints of tree rooting depth within a competitive environment, undermining the ecological realism of their results. Here, we present an approach of implementing variable rooting strategies and dynamic root growth into the LPJmL4.0 (Lund-Potsdam-Jena managed Land) dynamic global vegetation model (DGVM) and apply it to tropical and sub-tropical South America under contemporary climate conditions. We show how competing rooting strategies which underlie the trade-off between above- and below-ground carbon investment lead to more realistic simulation of intra-annual productivity and evapotranspiration and consequently of forest cover and spatial biomass distribution. We find that climate and soil depth determine a spatially heterogeneous pattern of mean rooting depth and below-ground biomass across the study region. Our findings support the hypothesis that the ability of evergreen trees to adjust their rooting systems to seasonally dry climates is crucial to explaining the current dominance, productivity and evapotranspiration of evergreen forests in tropical South America

Climate-induced hysteresis of the tropical forest in a fire-enabled Earth system model

Tropical rainforests are recognized as one of the terrestrial tipping elements which could have profound impacts on the global climate, once their vegetation has transitioned into savanna or grassland states. While several studies investigated the savannization of, e.g., the Amazon rainforest, few studies considered the influence of fire. Fire is expected to potentially shift the savanna-forest boundary and hence impact the dynamical equilibrium between these two possible vegetation states under changing climate. To investigate the climate-induced hysteresis in pan-tropical forests and the impact of fire under future climate conditions, we employed the Earth system model CM2Mc, which is biophysically coupled to the fire-enabled state-of-the-art dynamic global vegetation model LPJmL. We conducted several simulation experiments where atmospheric CO$$_2$$ concentrations increased (impact phase) and decreased from the new state (recovery phase), each with and without enabling wildfires. We find a hysteresis of the biomass and vegetation cover in tropical forest systems, with a strong regional heterogeneity. After biomass loss along increasing atmospheric CO$$_2$$ concentrations and accompanied mean surface temperature increase of about 4 $$^\circ$$C (impact phase), the system does not recover completely into its original state on its return path, even though atmospheric CO$$_2$$ concentrations return to their original state. While not detecting large-scale tipping points, our results show a climate-induced hysteresis in tropical forest and lagged responses in forest recovery after the climate has returned to its original state. Wildfires slightly widen the climate-induced hysteresis in tropical forests and lead to a lagged response in forest recovery by ca. 30 years

Earth beyond six of nine planetary boundaries

This planetary boundaries framework update finds that six of the nine boundaries are transgressed, suggesting that Earth is now well outside of the safe operating space for humanity. Ocean acidification is close to being breached, while aerosol loading regionally exceeds the boundary. Stratospheric ozone levels have slightly recovered. The transgression level has increased for all boundaries earlier identified as overstepped. As primary production drives Earth system biosphere functions, human appropriation of net primary production is proposed as a control variable for functional biosphere integrity. This boundary is also transgressed. Earth system modeling of different levels of the transgression of the climate and land system change boundaries illustrates that these anthropogenic impacts on Earth system must be considered in a systemic context